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Abstract. - Spatiotemporal control of excitable media is of paramount importance in the develop- 
ment of new applications, ranging from biology to physics. To this end we identify and describe a 
qualitative property of excitable media that enables us to generate a sequence of traveling pulses of 
any desired length, using a one-time initial stimulus. The wave trains are produced by a transient 
pacemaker generated by a one-time suitably tailored spatially localized finite amplitude stimulus, 
and belong to a family of fast pulse trains. A second family, of slow pulse trains, is also present. 
The latter are created through a clumping instability of a traveling wave state (in an excitable 
regime) and are inaccessible to single localized stimuli of the type we use. The results indicate that 
the presence of a large multiplicity of stable, accessible, multi-pulse states is a general property of 
simple models of excitable media. 



^-^ Excitable media are characterized by a large (finite am- 
QQplitude) response to a supra-threshold perturbation of 
^^the rest state, followed by decay back to the same rest 
^^tate [1]. Such behavior is frequently found in biologi- 
•^cal [2], chemical [3], and physical [4] systems. As sug- 
^S gested independently by FitzHugh and Nagumo [5] , tem- 
^poral excitable behavior can be understood qualitatively 
via a prototype two-component ordinary differential equa- 
tion of Bonhoeffer-van der Pol type with a cubic nonlin- 
earity in the activator field [1]. In spatially extended ex- 
citable media with activator diffusion, such a perturba- 
tion results in the formation of a solitary traveling wave, 
hereafter referred to as a traveling pulse [1]. Since the 
medium returns to the rest state after the passage of the 
pulse, repeated stimulation is required to generate a se- 
quence of pulses [6,7]. Here we identify a novel property 
of excitable systems of this type that allows us to generate 
economically sequences of traveling pulses through a one- 
time initial stimulus, whose shape controls the number of 
pulses within the wavetrain. 

In a comoving frame, the one-dimensional (ID) profile 
of a traveling pulse corresponds, in many systems, to a 
Shil'nikov homoclinic orbit [8-12]. Such orbits are known 
to be accompanied by (an infinite number of) additional 

(^) Present address: Department of Chemical Engineering, Tech- 
nion - Israel Institute of Technology, Haifa 32000, Israel. 



multi-pulse homoclinic orbits [7,9,13], and these in turn 
correspond to spatially localized groups of pulses (here- 
after finite pulse trains) traveling together with common 
speed. This speed differs in general from the speed of 
a single pulse. All of these states form close to the pa- 
rameter values at which the primary homoclinic orbit is 
present [9,11] and have similar speeds. In contrast, there 
are also spatially extended states (hereafter infinite pulse 
trains) consisting of copies of the single pulse state, in 
which each pulse is locked to the oscillatory tail of the 
preceding pulse [10,14]. Many different states consisting 
of equally or unequally spaced pulses are possible, each 
traveling at its own speed [6, 7, 10, 13, 15]. Although so- 
lutions of either type are readily constructed using ideas 
from spatial dynamics, their stability properties are in gen- 
eral unknown [7,10,12,13,15]. Moreover, even when the 
solutions are stable, no robust procedures are known for 
generating a pulse train of a desired type or length without 
continued input. 

In this Letter, we demonstrate that in spatially ex- 
tended homogeneous excitable media, distinct families of 
multi-pulse states can be simultaneously stable, and ex- 
plain how the different states may be generated. Of these 
the fast pulse trains can be generated with great selectivity 
by suitably tuned one-time perturbations, via the forma- 
tion of a transient pacemaker. We show how to mold the 
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initial perturbation to achieve the desired result. The ba- 
sic idea is simple: a one-time perturbation evolves first in 
amplitude and so forms a localized standing oscillation, 
i.e., a pacemaker. If this oscillation is in turn unstable 
to decay into solitary waves, the pacemaker is transient 
and the number of pulses generated by the perturbation 
is finite. In contrast, the slow pulse trains are the result 
of a long wavelength "clumping" instability of the travel- 
ing waves state. The results are obtained in one spatial 
dimension, but qualitatively similar results hold in two 
dimensions. 

Model Equations. — In two-component models an 
oscillatory instability of a homogeneous state necessarily 
comes in with a zero wavenumber at onset and hence pro- 
duces spatially uniform oscillations [16]. This is not the 
case for three-component models [17]. In the following we 
focus on a system of this type that is known to support, 
under identical conditions, both traveling solitary pulses 
and spontaneous pacing activity [18]: 

dtu = u-u^-v^D{dxx^dyy)u, 
dtv = Cy {u — ayV — au)W — ao) , 
dtw = ew{u-w)^ {dxx + dyy) w . 

Here u{x^y^t) can be thought of as an activator and 
v{x^ y, t) as an inhibitor; the latter is in turn controlled by 
a third component w{x^y^t)^ also generated in response 
to the activator. When a^j = eqs. (la) and (lb) reduce 
to the standard FitzHugh-Nagumo (FHN) form [1]. In 
what follows we set a,,, = 0.5, and adopt the parameters 
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Fig. 1: Bifurcation diagram for spatially homogeneous steady 
([/o,[/±, thin black lines) and oscillatory (heavy light colour 
lines) solutions of eqs. (1). Oscillatory branches of standing 
waves (SW) and traveling waves (TW) with fixed spatial period 
Ac = 2ti /kc — 4.07 are shown in terms of the maximum value 
Umax of u(x,t), and emerge from a finite wavenumber Hopf 
bifurcation of Uo = (uq^vq^wq) at ay = a^ — 0.4497. Solid 
(dashed) lines indicate stable (unstable) solutions. The shaded 
region, 0.4497 ^ a^ < 0.55, delimits the existence of stable 
pulse trains of fast type. The inset represents the dispersion 
relation at onset (a^ = a^ , kc — 1.54, ujc — 0.4, with TW 
phase velocity cj^ :^ 0.26) and near the rightmost limit of the 
shaded region, ay = 0.55, where Uo is linearly stable. 



shown in fig. 1, both branches bifurcate subcritically (i.e., 
w 



1.0, ao = -0.1, D = 0.005, with a^ used as in the a^ > a^ direction) but eventually turn towards 



a control parameter. 

Uniform and Nonuniform Solutions. — Equa- 
tions (1) admit three spatially homogeneous solutions re- 
ferred to as Uo , U- , [/+ ; the U± arise through a saddle- 
node bifurcation at a^ — a^^ :^ 0.8919 and are present 
for ay > a^^ only (see fig. 1). In excitable media the rest 
state, here Uq = {uo^vq^wq)^ has to be stable to infinites- 
imal perturbations. 





at-\-ikx 



C.C. 



o{s^ 



(2) 

Here a is the (complex) perturbation growth rate, k > 
is the wavenumber, and c.c. denotes a complex conjugate. 
Linear stability analysis (5 <C 1) shows that Uq is stable for 
: 0.4497 (the excitable regime) but unstable for 



> a, 



w 

V 



< a^ (see inset in fig. 1). The bifurcation at ay 



,w 



is a Hopf bifurcation with finite wavenumber kc — 1.54 at 
onset. The corresponding critical frequency is c^c — 0.4. 
Standard theory [19] now shows that two branches of oscil- 
latory states with spatial period Ac = 2ti /kc — 4.07 bifur- 
cate from the rest state at ay = a^, a branch of traveling 
waves (TW) and a branch of standing waves (SW) ^. As 



small ay] for SW, the subcritical regime is very small. It 
follows that near onset the TW are once unstable, while 
the SW are twice unstable. Moreover, the TW acquire sta- 
bility with respect to wavelength Ac perturbations once 
the TW branch turns around towards smaller a^, while 
the SW remain unstable to TW perturbations. We shall 
see, however, that the TW are in fact unstable beyond the 
saddle-node {ay :^ 0.63) with respect to long wavelength 
perturbations (fig. 1). These instabilities accumulate at 
ay ^ 0.5 so that the TW are stable with respect to all 
perturbations for ay < 0.5 only, cf. [21]. 

Generation of Pulse Trains. — In addition to the 
spatially periodic SW and TW states eqs. (1) also admit 
different types of single-pulse and multi-pulse solutions. 
In a comoving reference frame these states correspond to 
homoclinic orbits, i.e., to solutions asymptotic to Uq at 
^ ^ ±oo. There are in fact two families of such states 
referred to in what follows as slow and fast, based on their 



frame the TW corresponds to a periodic orbit with period Ac, with 



the speed c : 



of the waves computed as a nonhnear eigenvalue 



-"^The TW were computed in the comoving frame ^ = x - 



i.e., by replacing dt in eqs. (1) by 



d^, and dx by d^. In this 



using the numerical continuation package AUTO [20]. Temporal 
stability was computed via a standard numerical eigenvalue method 
using the time-dependent version of eqs. (1) in the comoving frame. 
The SW branch was computed via direct integration of eqs. (1) in < 
X < Ac/4 with Neumann boundary conditions at x = and Dirichlet 
boundary conditions at x = Ac/4, to suppress the instability of the 
SW with respect to TW perturbations. 



Generation of finite wave trains in excitable media 




-aII|L: 2x^ 




--iW : 3\ 




L_ 


\i \J \l Vl 



4X 



Fig. 2: The "clumping" instability of traveling waves and 
the formation of finite pulse trains at av = 0.6 in small pe- 
riodic domains with period (a) L = 2 Ac — 8.14, (b) L = 
3Ac ^ 12.21, and (c) L = 4Xc ^ 16.28. The background 
state filling the domain corresponds to the homogeneous state 
Uo = {uo,vo,wo) c^ (-0.650,-0.375,-0.650). Left panels: 
space-time plots with x along the vertical axis and time hor- 
izontally (dark colour indicates larger values of the u field). 
Right panels: profiles of finite pulse trains at t = 390, travel- 
ing without change towards the right. 



speed of propagation relative to that of the TW. 

Slow pulse trains. Formation of the slow pulse trains 
is associated with the long wavelength instability of the 
TW (fig. 1). These instabilities are present in the regime 
0-5 ^ ay < 0.63 and the number of slow pulses that re- 
sults depends on the length L of the periodic domain used 
to solve eqs. (1). For example, for a^ = 0.6 a bound 
state of two pulses is generated when the domain period 
is 1/ = 2 Ac, while groups of three and four pulses form 
when L = 3 Ac, and L = 4Ac, respectively (fig. 2). These 
pulse trains propagate with speeds c o^ 0.08 that are sub- 
stantially slower than c^^ c^ 0.17, as indicated by the 
dramatic change in slope of the space-time trajectories of 
each crest as the instability evolves. The corresponding 
solution profiles are shown in the right set of panels in 
fig. 2. No single traveling pulses are generated for this 
parameter choice, presumably because the TW are stable 
in an 1/ = Ac periodic domain. 

We may think of the resulting instability as a "clump- 
ing" instability, and believe that it is a consequence of the 
fact that the fastest growing mode in an L = nAc domain 
is the mode with wavelength nAc- This property is char- 
acteristic of a phase-like instability called the Eckhaus- 
Benjamin-Feir instability except that here this instability 
does not result in phase slips [22]: instead when the wave 
amplitude falls locally to sub-threshold value the system 
nucleates the rest state Uq. The result is evolution into a 
clumped state that conserves the number of crests. 

Fast pulse trains. The fast pulse trains, on the other 
hand, are generated not by linear instability but in re- 
sponse to a one-time finite amplitude spatially localized 
perturbation. Such perturbations evolve more rapidly in 
amplitude than in phase and hence result in a spatially lo- 



calized standing oscillation. These, like the spatially peri- 
odic SW, are in turn unstable to TW perturbations. Thus 
in the shaded region in fig. 1, suitably chosen finite ampli- 
tude perturbations evolve into transient pacemakers whose 
decay time controls the length of the pulse train that is 
emitted. As a result the number of pulses in the wavetrain 
is now controlled by the decay time of the pacemaker and 
hence the shape of the initial perturbation, instead of the 
imposed spatial period, L. 

To demonstrate the development of distinct but finite 
pulse trains in large domains {L ^ Ac) by the above mech- 
anism, we focus on the parameter regime in which both 
Uo and TW are linearly stable, i.e., 0.45 ^ ay < 0.5. In 
what follows, we consider eqs. (1) with Neumann bound- 
ary conditions on the ID domain < x < L, and examine 
the effect of different initial perturbations on pacemaker 
emergence and the resulting solutions. The perturbations 
have a smoothed top-hat profile centered at x = 0, 



u{x,t = 0) = uo^—(l 



tanh - 



R 



(3) 



where A denotes the amplitude of the perturbation, R <^ 
L its width, and A < R measures the steepness of the 
initial interface. When ay > a^ and A ~ 0.1, i? ~ 1, the 
excitability threshold is locally exceeded and the system 
responds with a single propagating pulse. Figure 3(a), 
shows the propagation of such a pulse in a space-time plot 
(left panel) with an enlargement of the initial phase (mid- 
dle panel) and the spatial pulse profile (right panel). The 
figure shows that the initial perturbation grows abruptly 
in place before decaying into a propagating pulse. Once in 
motion, the profile of the pulse becomes distinctly asym- 
metrical, with a train of spatially decaying oscillations to- 
wards its rear, as expected of a Shil'nikov saddle- focus [8]. 

While such single pulse profiles are known from the 
two- variable FHN system [1,12], the three- variable model 
[eqs. (1)] allows the generation of more complicated solu- 
tions such as A/"- homo clinic orbits of Shil'nikov type [9]. 
To observe these states we decrease the stimulus width 
to R = 0.4. This initial perturbation sets up a tran- 
sient pacemaker and the system now generates a finite 
pulse train [fig. 3(b), left panel], even though no addi- 
tional stimulation is applied after t = 0. The middle panel 
of fig. 3(b) shows that the pacemaker oscillates in the form 
of a standing oscillation in < x < 0.25 for three cycles, 
before decaying. Each oscillation cycle that exceeds the 
excitability threshold results in the generation of a travel- 
ing pulse. Consequently an initial stimulus with A = 0.3, 
R = 0.4 triggers three pulses moving together as a group 
[see fig. 3(b)]. 

Numerically, we find that the single pulse in fig. 3(a) 
propagates with speed c^^^*^ :^ 0.2482, while the pulse 
trains all propagate with essentially identical and slightly 
larger speed c :^ 0.2487, to within numerical error of order 
of 10~^. Both velocities are substantially larger than the 
phase velocity of the TW, c^^ c^ 0.1877 at ay = 0.5. 
Consequently we call the resulting finite pulse trains fast. 
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Fig. 3: Solutions of eqs. (1) at a^ = 0.5 with Neumann boundary conditions at x = 0, x = 40 using different stimuli at t = 0: (a) 
A — 0.3, R — 1^ (b) A — 0.3, R — 0.4, both for A = 0.1. Left panels: space-time plots with x along the vertical axis and time 
horizontally (dark colour indicates larger values of the u field). Middle panels: enlargement of the space-time plots in vicinity 
of the initial perturbation, x G [0, 1], demonstrating (a) the initial growth of the disturbance in amplitude from u = uq + A <^ 
to 1^ ~ 1, followed by its decay when the initial perturbation is broad, and (b) the formation and eventual decay of a transient 
pacemaker when the initial perturbation is narrower; the dashed line in (b) is at x = 0.25. Each oscillation of the pacemaker 
generates one pulse. Right panels: profiles of (a) single-pulse and (b) three-pulse states, both traveling without change towards 
the right. 



Comparison with fig. 2 shows that the fast pulses are much 
broader than the slow pulses, but otherwise resemble each 
other. 

The finite pulse trains accessible by varying the ampli- 
tude A and width R of the initial perturbation are summa- 
rized in fig. 4. In the vicinity of the extinction boundary 
the number of pulses increases, albeit in a nonmonotonic 
fashion, and this property of the system can be exploited 
to identify empirically the initial conditions that result 
in the desired number of pulses. In the bottom panel of 
fig. 4 we provide a more detailed sample of the possible 
outcomes. Evidently the number of pulses is very sensitive 
to the form of the initial perturbation, and this sensitiv- 
ity extends to the dependence of the results on the length 
scale A introduced in eq. (3). 

Qualitatively similar results are found for other pa- 
rameterized families of initial conditions, provided these 
broadly resemble the top-hat profile in eq. (3). 

Bistability between Different Pulse Trains. — 

The similarity in the pulse profiles suggests that the 
two families of finite pulse trains may in fact be related. 
This supposition is supported by our numerical simula- 
tions, which show that stable pulse trains of slow and 
fast type may coexist, cf. [18]. To identify the connec- 
tion between these two families we have again employed 
numerical continuation [20] but this time initiated with 
single, double and triple homoclinic profiles of fast type, 
such as those shown in fig. 3. In fig. 5, we show that the 
slow and fast pulse trains are in fact separated by un- 
stable branches (dashed lines). This result explains why 
we find two families of traveling pulse trains, and why we 
observe bistability between them. Observe that the fast 
pulse trains all travel with essentially the same speed c as 
ay decreases, and do so over a broad interval of a^ val- 
ues. This is the case for all the multi-pulse states shown 
in fig. 4. This fact facilitates the insertion of additional 



R 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 






























- 





























- 








3 


3 \ 


'3" 


3 








1 





3 


3 


3 ; 


2 










1 - 





3 


3 


3 


2 










1 - 


- 3 


3 


3 


3 


2 










1 - 


3 


3 


3 


3 


2 










1 - 


3 


3 


3 


2 


2 










1 - 


3 


3 


3 


2 


2 










1 - 


- 3 


3 


3 


2 


2 










1 - 





0.46 


0.48 


0.50 


0.52 


0.54 


R 
0.56 


0.58 


0.60 


0.62 


0.64 0.66 


0.3 


3 


3 


3 


3 


3 


3 


3 


3 


6 


4 1 - 


0.32 


- 3 


3 


3 


3 


3 


6 


4 


6 


4 


1 1 - 


0.34 


3 


3 


3 


3 


4 


2 


2 


9 


4 


1 1 - 


0.36 


3 


3 


3 


4 


2 


2 


5 


8 


1 


1 1 - 


0.38 


3 


3 


4 


2 


2 


2 


6 


4 


1 


1 1 - 


0.4 


3 


5 


2 


2 


2 


2 


9 


1 


1 


1 1 - 



Fig. 4: Top panel: finite fast pulse trains computed from the 
initial perturbation (3) with A = 0.1 for different values of A 
and R when a^ = 0.5. The integers correspond to the number 
N of pulses in the wavetrain, the bold italic fonts correspond- 
ing to the solutions described in fig. 3. Bottom panel: variation 
in N in the vicinity of the extinction boundary, in the region 
indicated by the small rectangle in the top panel. The numeri- 
cal integration of eqs. (1) was performed with Ax = 0.005 and 
Neumann boundary conditions on an L = 50 domain. 



pulses into a pulse train through incremental changes in 
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Fig. 5: The speed c of the fast (heavy colour lines) and 
(thin colour lines) pulse trains as a function of ay tog( 
with the phase speed of traveling wave (TW) states. Only 
three branches of traveling pulse trains are shown. The phase 
speed of the TW lies between the speeds of the fast and slow 
pulse trains. Solid (dashed) lines indicate stable (unstable) 
states. The accumulation of the saddle-nodes of the fast pulse 
trains (a^ ~ 0.55) defines the upper boundary of the shaded 
region displayed in fig. 1. The pulse trains were computed in 
a comoving frame, on a periodic domain with L — 100. 



the initial stimulus. The resulting selectivity is a conse- 
quence of the large multiplicity of coexisting stable pulse 
trains and the intertwined structure of their basins of at- 
traction. In contrast, the branches of slow pulse trains all 
turn around towards smaller a^ as c decreases (not shown), 
and in so doing follow the behavior of the TW branch. 

The results in fig. 5 help us to understand the observed 
absence of fast pulse trains at a^ = 0.6, and suggest that 
the empirically determined upper boundary of the shaded 
region in fig. 1, i.e., the region within which finite am- 
plitude perturbations generate finite pulse trains, is de- 
termined by the accumulation point of the saddle-node 
bifurcations at which the fast pulse trains lose stability 
and cease to exist. As seen in the figure these saddle- 
node bifurcations accumulate at ay c^ 0.55, an observa- 
tion that can be confirmed by numerical continuation of 
the additional multi-pulse states shown in fig. 4. This 
value is close to that at which the TW solutions acquire 
stability with respect to wavelength nXc (n > 1) perturba- 
tions but is not related to it. Numerically we find that in 
(^v ^ 0-55 the finite amplitude perturbations always reach 
fast pulse trains, even when these coexist with stable slow 
pulse trains. Figure 5 also shows that for ay > 0.55 slow 
pulse trains, including the missing slow single-pulse state 
in L = Ac domains, should be reachable by appropriate 
finite amplitude perturbations of Uq. 

Generation of Finite Trains of Target Waves. — 

We now turn to the corresponding results in two spatial 
dimensions (2D), with Neumann boundary conditions on 
a square domain. The initial perturbation is taken to be 
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Fig. 6: Solutions of eqs. (1) in the {x,y) plane with ay = 0.5, 
using different initial stimuli. Dark colour implies larger values 
the u field; the propagation direction is from the top- left corner 
to the bottom-right corner, (a) Single wave propagation for 
R = 2.0. (b) Transient pacing and the formation of a finite 
number of target waves for R — 0.8. (c) Formation of a stable 
pacemaker and persistent target wave emission for i? = 0.6. In 
all cases the perturbation amplitude is A = 0.3, with A = 0.1; 
the physical domain size is 20 x 20 with Neumann boundary 
conditions, and Ax = Ay — 0.005. 



spot-like and centered at the top-left corner of the do- 
main, as shown in fig. 6; the parameter R now measures 
its extent in the radial direction. 

The two types of behavior, propagation of a single tar- 
get wave [see fig. 6(a)] and of a finite train of target waves 
due to transient pacing activity [see fig. 6(b)], persist also 
in 2D; the spatial cross- sections of these states along the 
propagation diagonal [figs. 6(a,b)] admit profiles that re- 
semble the profiles of fast family type obtained in ID (see 
fig. 3). The value of R required for the finite target waves is 
now larger {R ^ 1) than in the ID case, but with increas- 
ing R single-pulse propagation results, as in ID. However, 
in contrast to ID a one-time perturbation with a small 
value of R produces persistent pacing activity that results 
in standard target waves [see fig. 6(c)] [18]. These results 
are not affected by the choice of partial domain and apply 
equally to a large square domain with a spot-like pertur- 
bation in the center. We attribute the differences between 
the ID and 2D results to curvature effects in the pace- 
maker region. 

Conclusions. — We have demonstrated that stable 
multi-pulse solutions can be generated in a spatially ex- 
tended excitable medium by a one-time stimulus, with no 
need for additional stimulation. These finite pulse trains 
are the result of the formation of a transient pacemaker 
with decaying amplitude that generates solitary waves 
while its amplitude remains above threshold for pulse gen- 
eration. The primary requirement for this property of the 
system is a nonlinear (subcritical) instability of the uni- 
form state to both traveling waves (TW) and standing 
waves (SW), as in fig. 1. We have seen that under these 
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conditions stable multi-pulse states are present, and that 
the basins of attraction of these states can be selected by 
a suitably tailored initial stimulus through the excitation 
of a transient pacemaker. We believe these results to be 
independent of the details of the system used here, and 
thus characteristic of a broad class of excitable systems. 

In the present example the pulse trains generated by the 
transient pacemaker propagate with a speed substantially 
larger than the coexisting spatially periodic TW, and may 
also coexist with a distinct family of traveling pulses that 
propagate more slowly than the TW. The latter form as 
a result of a long wavelength "clumping" instability of 
the TW that is also a consequence of the subcriticality 
of the TW branch. However, the slow pulse trains are 
not accessible by finite amplitude perturbations of the rest 
state of the type we use. 

The study was conducted within the phenomenologi- 
cal framework of a three- variable FitzHugh-Nagumo type 
model [18], and employed bifurcation analysis, numerical 
continuation, and parametrized stimuli to generate a va- 
riety of finite pulse trains, as summarized in fig. 4. These 
results demonstrate that near the extinction boundary the 
length N of the fast pulse train becomes a sensitive func- 
tion of the shape of the initial stimulus, enabling us to 
"dial in" different pulse trains by appropriately shaping 
the initial stimulus. The details of the resulting look-up ta- 
ble (fig. 4) are system-specific but need only be computed 
once; similar look-up tables can be established empirically 
for each model of interest. As shown in figs. 6(a,b) these 
results extend to two space dimensions as well. 

We speculate that this intriguing property of excitable 
systems may provide new insights into a number of ap- 
plied science processes exhibiting reaction-diffusion behav- 
ior [1,17,23], including chemical reactions [3], excitation 
and propagation of signals in neurons and cardiac tissue, 
and intracellular calcium cycling [2]. Technological ap- 
plications to information transfer using nonlinear optical 
systems [24] and to unconventional computation based on 
reaction-diffusion systems [25] are also possible. 
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